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The quantum Heisenberg model belongs to the most intensively studied models in statistical mechanics 
(see [1] and references therein), since it exhibits very interesting quantum phenomena. However, due to the 
mathematical complexities related to the quantum nature of the model, only few exact results have been 
obtained in this field so far. Owing to this fact, the main aim of this work is to introduce a new statisti- 
cal model, which is exactly solvable for planar lattices and simultaneously exhibits interesting quantum 
properties. The model to be investigated here consists of Ising- and Heisenberg-type atoms regularly dis- 
tributed on planar doubly decorated lattices, as it is depicted in Fig. 1 for the case of the square lattice. In 
order to derive the exact results for the model under investigation, it is very convenient to express the total 
Hamiltonian in the form Tid — J2^2{ 2 ^k, where the summation is carried out over all the bonds of the 
original (undecorated) lattice, N denotes the total number of the Ising-type atoms and q is the coordination 
number of the original (undecorated) lattice. The bond Hamiltonian Hk is given by (see Fig. 1) 

Hk = -J ^{S% 1 S% 2 + ^Slz) + S^Slz - Ji(§ki + 5*fc2)(Afei + AD - D (5"fci) 2 + (^2) 

where SI (7 = x, y, z; a — 1, 2) represent the relevant components of the standard spin-1 operators and 
Mfca those of the spin- 1/2 operators. The exchange interaction couples the nearest-neighboring Ising 
and Heisenberg atoms and similarly, J couples nearest-neighboring Heisenberg pairs. Finally, A describes 
the anisotropy in the Heisenberg interactions and D denotes the single-ion anisotropy. Taking into account 




(1) 



Fig. 1 Part of the doubly decorated square lattice. Black circles represent the Ising atoms with spin 1/2 and gray ones 
denote the decorating Heisenberg atoms with spin 1 . 
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the standard commutation relations for the bond Hamiltonians (i.e. [Hi, Wfc] = 0, i ^ k), one can rewrite 
the partition function of the decorated model in the form 

Nq/2 

Z d = Trexp(-(3H d ) = Tr {jt(} Tr S fe i Tr s fc2 exp(-0H k ). (2) 

fe=i 

Here, the symbol Trr^-j. means a trace over all degrees of freedom of Ising spins, Tis kl TTs k2 denotes a 
trace over a couple of Heisenberg spins residing on the fc-th bond and = l/(k-QT). To proceed further, 
we introduce an extended decoration-iteration transformation [2] . Namely, 

Tr Sfcl Tr S)c2 exp(-/3?4) = Aexp{PRfi z kl fi z k2 ), (3) 

from which, after evaluating the l.h.s., we obtain for the transformation parameters A and R the following 
expressions: 



where 



A = y/(W + Wi)(W + W 2 ), 0R = 2 ln(W + Wi) - 2 ln(W + W 2 ), (4) 



Wo = exp[/?(2L> - J)] + 2 exp[/?(L> - J/2)] cosh 



2 V(2^D- J) 2 +8(JA)2 



(5) 



Wi = 2exp[/3(2/J + J)]cosh(2/?Ji) + 4exp(/3D)cosh(/3JA)cosh(/3Ji), (6) 
W 2 = 2exp[/3(2/J + J)]+4exp(/3/J))cosh( / 9JA). (7) 

Finally, after substituting (3) into (2) one obtains the equation Z d = A Nq / 2 Z , which relates the parti- 
tion function of doubly decorated Ising-Heisenberg model (Z d ) to that one of the original (undecorated) 
spin- 1/2 Ising model (Zq), that is exactly known for many planar lattices. Although, from this simple re- 
lation we can directly calculate all relevant physical quantities over the whole range of the parameters and 
temperature, we restrict ourselves in this brief report only to the analysis of the ground-state properties of 
the system. For this purpose, we have studied the lowtemperature behavior of the sublattice magnetization 
rriA = (/*fe a )i m b = (Ska)' as we ^ as trie variations of the relevant nearest-neighbor correlation func- 
tions, namely, qff = (fi z kl fi z kl ), q% = (Ali^fi). qll = (S^S^) (7 = x, y, z) and 7] = {(S z kl ) 2 ). Before 
discussing the results, it is worth noticing that the ground state of our system is universal in the sense that 
it does not depend on the coordination number and the space dimensionality of the lattice. 

At first, we have depicted in Fig. 2a the ground-state phase diagram in the D — Jl space for the 
system without single-ion anisotropy, i.e. D/J = 0.0. As one can see, three different phases can appear 
in the system depending on the ratio between J\ and D. The phase denoted as Oi represents a standard 
ferromagnetic phase which is usually observed in the pure Ising systems. As one can expect, in this case 
we have m A = 0.5, m B = 1.0, q zz = 0.25, q z f- = 0.5, q z h z h = 1.0, rj = 1.0 and q%% = qf h = 0.0. On the 
other hand, the phase denoted as O2 exhibits some unusual features, though it is also a long-range ordered 
ferromagnetic phase in which all Ising spins are perfectly aligned along the easy axis (thus, tua — 0.5 and 
iti = 0.25). Contrary to the perfectly ordered Ising sublattice, the Heisenberg one behaves in a different 
way. Indeed one finds that tub — 0.5, q\i = 0.25, q 7 ^ h — 0.0 and r\ — 0.5. Moreover, in this phase one also 
observes the nonvanishing short-range order in the xy plane (qf^ = qW — 0.5). This rather unexpected 
behavior is simply explainable within the valence-bond-solid picture [3], that is schematically depicted for 
the decorating atoms in the inset of Fig. 2a. This is, of course, the pure quantum effect which is caused by 
relatively strong quantum fluctuations in the system. It is also clear that the phase O2 represents a mixture 
of regularly ordered Ising spins with the quantum phase on the Heisenberg sublattice and to the best of our 
knowledge such a phase has not been discussed in the literature before. Finally, if the anisotropy parameter 
A becomes strong enough, a disordered phase (DP) appears in the system. In addition to the trivial fact 
that in this phase mA = tub = 0.0 is found, one also finds q zz = 0.25, q z t z — q z t z — 0.0, q^ z h < 0.0 and 
Qhh ~ Ihji > 0-0. In fact, all these findings clearly indicate that the relevant disordered quantum phase is 
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Fig. 2 Ground-state phase diagrams of the doubly decorated Ising-Heisenberg model. Oi means the standard ferro- 
magnetic phase, O2 represents the quantum ferromagnetic phase, the properties of which can be explained within the 
valence-bond-solid picture. The situation is schematically depicted for the decorating atoms in the inset of Fig. 2a. 
The large open circles represent the atomic sites with S = 1, the small solid circles and arrows denote the spin 1/2 
variables. The solid line shows the valence bond. Finally, DP denotes the disordered frustrated phase and the black 
points mark in Fig. 2b the triple points. 



strongly frustrated. It is also worth noticing that the boundary between the phase Oi and O2 is given by 
Ji = J(A - 1) and that separating the phase 2 and DP is given by J x = J(0.5\/l + 8A 2 - A - 0.5). 

To complete the ground-state analysis, we have shown in Fig. 2b the phase diagram of the model in 
the A — D space for J\ — J. Similarly, as in the previous case we have again observed the same phases. 
The phase diagram, however, becomes more interesting, since in this case we have also found two triple 
points, namely, Ti = [2 — \/2, — s/% and T 2 = [2 + \/2, \/2] . Moreover, it is apparent that the relevant phase 
boundaries become more complex. Here we only state that in both cases (D = and D ^ 0) all phase 
boundaries represent the points at which the system exhibits first-order phase transitions. The detailed 
analysis of the ground state properties, as well as that one of the thermal behavior of the system is being 
left for our future work which is now in preparation. 

To conclude this work, we would like to emphasize that the Ising-Heisenberg models (including the 
present model) represent very useful exactly solvable statistical models that can exhibit very interesting 
quantum phenomena. Moreover, these models can be helpful in description of some real materials. For 
example, the recently synthesized compound studied in Ref. [4] has the structure of the doubly decorated 
honeycomb lattice. 
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